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Abstract 

We present the holographic predictions for cosmological 3-point correlators, involving 
both scalar and tensor modes, for a universe which started in a non-geometric holographic 
phase. Holographic formulae relate the cosmological 3-point functions to stress tensor 
correlation functions of a holographically dual three-dimensional non-gravitational QFT. 
We compute these correlators at 1-loop order for a theory containing massless scalars, 
fermions and gauge fields, and present an extensive analysis of the constraints due to Ward 
identities showing that they uniquely determine the correlators up to a few constants. We 
define shapes for all cosmological bispectra and compare the holographic shapes to the 
slow-roll ones, finding that some are distinguishable while others, perhaps surprisingly, are 
not. 
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1 Introduction 



In a recent series of papers [H El El HI El E] we put forward a holographic framework for infla- 
tionary cosmology and discussed a novel class of models describing a universe that started in a 
non-geometric phase, which is described holographically via a large- three-dimensional QFT. 
The power spectra and the scalar bispectrum were computed in [H |2] and [3], respectively, 
and in this paper we complete this program by computing the non-Gaussianities that involve 
tensors. Non-Gaussianities involving tensors are not expected to be measurable in near-future 
experiments. Nevertheless, they are still interesting theoretically and their structure has been 
the topic of several recent papers [3 El El ED] ■ 

The holographic model is specified by providing the dual QFT and the holographic dic- 
tionary that relates QFT correlation functions to cosmological observables. We worked out 
the holographic dictionary for non-Gaussianities involving tensors in [6], and in this paper 
we compute the relevant QFT correlation functions. The models we discuss are based on 
perturbative three-dimensional QFTs that admit a large- limit and have a generalised con- 
formal structure |lHll2j. An example of such a theory is SU {N) Yang-Mills theory coupled to 
massless scalars and fermions, with all fields transforming in the adjoint of SU{N). The non- 
Gaussianities are extracted from the 3-point function of the stress tensor of this theory. The 
leading order 1-loop computation of this 3-point function is independent of the interactions of 
the QFT, and thus our main task is to compute this 3-point function for free QFTs. 

Since all leading order results depend only on the free theory, let us briefly discuss the case 
in which the holographic model is a free QFT. In such a model, the spectrum is the exactly 
scale-invariant Harrison-Zel'dovich spectrum and the bispectrum is given exactly by the results 
reported here, i.e., the leading order results are the exact answer in the free theory. The shapes 
associated with the bispectrum may thus be considered as the analogue of the exact scale- 
invariant spectrum for higher point functions. We have seen in ^ that the scalar bispectrum 
shape for this model is indeed special: it is exactly equal to the factorisable equilateral shape, 
originally introduced in [13j. One may thus anticipate that shapes associated with the other 3- 
point functions will also have special properties. Possible shapes for the bispectrum involving 
only tensors have been discussed recently in [7] and here we will define and discuss shapes for 
the bispectrum involving both tensors and scalars. 

The computation of the 3-point function of the stress tensor at 1-loop is a non-trivial 



^ The holographic model is also the only model that yields exactly this shape. 
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task even for free QFTs. We discuss and develop several methods for evaluating the relevant 
Feynman diagrams. The 1-loop result is constrained by Ward identities and these provide a 
very non-trivial check of the expression we obtained by a direct computation. 

As mentioned above, to 1-loop order only the free part of the QFT enters. The QFT 
consists of gauge fields, fermions, minimal and conformal scalars. Conformal scalars and 
fermions are conformal field theories and their 3-point functions are constrained by conformal 
Ward identities. As is well known (from a position space analysis) [13], the 3-point function of 
the stress tensor in d = 3 is uniquely fixed by conformal invariance to be a linear combination 
of two conformal invariants, and is thus parametrised by two constants. (We assume parity 
is preserved). Our computation is done in momentum space and we thus provide the most 
general such 3-point functions in momentum space, where the two parameters are the number 
of conformal scalars and the number of fermions. The same computation was also recently 
reported in [t]^. We have explicitly verified that our results satisfy the conformal Ward 
identities. 

Let us now turn to minimal scalars and gauge fields. In three dimensions vectors are dual 
to scalars, so one may expect that gauge fields contribute the same as minimal scalars at 
1-loop order. We will indeed verify that this is the case. Note that beyond 1-loop the two are 
expected to contribute differently. Minimal scalars differ from conformal scalars in the way 
they couple to gravity, which in flat spacetime is reflected in their having a different stress 
tensor. More precisely, the stress tensor Tf- for a minimal scalar may be decomposed into a 
part Tf- corresponding to the stress tensor for a conformal scalar plus an "improvement term" : 

4=^-\{5^,^'-^.^,)<|?. (1.1) 

It follows that the 3-point function of Tf- may be computed from the 3-point functions involv- 
ing Tf- and the dimension one operator Oi = </>^. In turn, these 3-point functions are uniquely 
determined by conformal invariance, up to constants [14 1. Thus, effectively all 3-point func- 
tions are determined by conformal 3-point functions at this order, even though the underlying 
theory is not conformal. 

In [7J, the 3-point functions for tensors were computed in a de Sitter background. The 
de Sitter isometries act as the conformal group at late times, and the 3-point functions are 
then constrained by conformal invariance to be speciflc linear combinations of the 3-point 
functions of conformal scalars and fermions. Note that the de Sitter result is the leading order 
^Our results agree with the ones in v2 of 0. Relative to [7], we also computed semi-local terms. 
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approximation for slow-roll inflation. In general one expects that (broken) conformal invari- 
ance would constrain cosmological correlators in asymptotically de-Sitter slow-roll inflation, 
see also [151 US] • 

Our holographic results are for a very different universe, but we have seen that all relevant 
3-point functions are essentially determined by conformal 3-point functions. One may then 
wonder how our results compare with those of slow-roll inflation. The 3-point functions involv- 
ing only tensors are determined by the 3-point functions of conformal scalars and fermions, 
and, as in the discussion of [7], they agree exactly with slow-roll inflation if the fleld content 
of the dual QFT is appropriately chosen. The other 3-point functions (involving both scalar 
and tensor perturbations) are different but, perhaps surprisingly, they are rather similar. To 
quantify the differences we define (and plot) shape functions for all correlators, generalising 
the notion of shape functions for 3-point functions of only scalars or only tensors. 

This paper is organised as follows. In Section [2] we discuss the dual QFT and in Section [3] 
we present the holographic dictionary. In Section [J] we compute all relevant QFT correlation 
function by direct evaluation of the relevant Feynman integrals, and in Section [5] we explain 
their structure using Ward identities. The holographic predictions for the cosmological ob- 
servables are presented in Section O and these results are compared with the slow-roll ones in 
Section [71 We discuss our results in Section [51 Several technical results are presented in four 
appendices: in Appendix IX\ we summarise our notation and conventions for the helicity ten- 
sors, in Appendix IbI we present three different methods for evaluating the relevant diagrams, 
in Appendix [O we show that ghosts and gauge fixing terms do not contribute in correlators 
of the stress tensor and in Appendix [Dl we present the conformal and diffeomorphicm Ward 
identities. 



2 Dual QFT 

As a dual QFT we consider super-renormalisable theories that admit a large N limit and con- 
tain one dimensionful coupling constant. A prototype examplqj is three-dimensional SU{N) 
Yang-Mills theorj] coupled to a number of massless scalars and massless fermions, all trans- 
forming in the adjoint of SU{N). Theories of this type are typical in AdS/CFT where they 
appear as the world volume theories of D-branes. A general such model that admits a large N 

different exampie would be to consider 0{N) models, see ^7] for a related discussion. 
*We use the unconventional notation SU{N) as we reserve A'' for the analytically continued value iV — —iN. 
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limit is 

S = ^ [ d^x tr (jF.^jFIj + Ud(l)-^ f + Udx^f + i^^0i^^ + interactions ) , (2.1) 

where for all fields, = ip^'T'', and trTT^ = (5"^ We work with the Wick rotated QFT (of sig- 
nature (+, +, +)). This is mostly for convenience; we could equally well have stated all results 
in Lorentzian signature. The analytic continuation relevant for cosmology (which will appear 
in (|3.8p below) is a different continuation: it acts on the magnitude of the momentum. The 
gamma matrices satisfy {7j,7j} = —26ij. We consider Ma gauge fields {I = 1, . . . , Ma), 
M(f, minimal scalars 0'^ (J = 1, . . . , M^), M^ conformal scalars = 1) • • • ) -^x) ^^"^ 

M^ fermions (L = 1, . . . , M^). Note that Qym '^^^ dimension one in three dimensions. 
In general, the Lagrangian (|2.ip will also contain dimension- four interaction terms (see [2J). 
We will leave these interactions unspecified, however, as they do not contribute to the leading 
order calculations we perform here. 

In the next section we will present the holographic formulae that relate cosmological 3-point 
functions to correlation functions of the dual QFT. Generally speaking, the terms appearing 
in these formulae (see (I3.3p - (l3.6p ) are either 3-point functions of the stress tensor, or else semi- 
local terms (i.e., terms which are analytic in two of the three momentum). The semi- local 
terms involve either 2-point functions of the stress tensor or the T tensor defined by coupling 
the QFT to gravity, differentiating the stress tensor w.r.t. to the background metric and then 
setting the background metric to the flat metric. 



5Tij{xi) 



5g^\x2 



6^S 



ijki{xi,X2) = ^Jif^-^ „ = "^YJTT^T^UdT^ Q + 2^ij{^i)^kiKxi - ^2)- (2.2) 



5g'^^ {xi)8g^'-{x2 



1. 



The 2-point function takes the following general form 

{{Tij{q)TM{-q))) = A{q)U^Jkl + B{q)7rijTTM, (2.3) 

where the double bracket notation suppresses the delta function associated with momentum 
conservation, i.e., 

m,{^,)TM{q2)) = {27TfS{^i + q2){{Tij{qi)ni{-qi))), (2.4) 
and the transverse and transverse traceless projection operators are respectively 

TTij = 5ij - , 'n.ijki = -^^{'n'ik'^ji + -n-wn-jk - T^ijT^ki)- (2.5) 
q z 
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T(qJ 



Figure 1: 1-loop contribution to 
the stress tensor 3-point function. 



The leading contribution to the 2- and 3-point function comes from 1-loop diagrams (see 
Fig. [T]for the 3-point function), which are of order iV^ and involve only the free part of the 
Lagrangian. Interactions contribute to diagrams at 2-loop order and higher, but these are 
suppressed by factors of g^g relative to the 1-loop contribution and will be neglected here. As 
discussed in [Il[2], g^s is of the order of n^-l ~ O(10~^). Indeed, fitting the WMAP data to 
this model [5] leads to a small value of g'^g justifying the perturbative treatment. 

For spatially flat cosmologies, the background metric seen by the dual QFT is also flat. 
The dual stress tensor is then given by 

o xq 

T-.-——— - A. T'^ A. T"^ A. T^- (2 6) 

~ x/gSg'^ ~ 

V 9ij—°ij 

where the contributions from the various fields (suppressing the interactions, as well as the 
ghost and gauge-fixing terms which we discuss in Appendix E]) in (|2.ip are 

1 . r„r„r .1 



rpX 



2 " [-^i ^ ik - ^ij jFtiFi 



kl-^kl] ; 



1 



9ym 



tr[di 



1 



9yu ° 2 8 



1 



■ 1 



9ym 



(2.7) 
(2.8) 
(2.9) 
(2.10) 



where d 



di — di and symmetrisation is performed with unit weight. Note that the trace 
of the stress tensors for both conformally coupled scalars and for massless fermions vanish on 
shell. This is a consequence of the Weyl invariance of the quadratic action for these fields 
(with the fields transforming non-trivially) when the action (j2.ip is appropriately coupled to 
gravity. 



6 



3 Holographic formulae for cosmological 3-point functions 



We review in this section the holographic formulae for cosmological 3-point functions derived 
in [6]. These formulae relate the late-time behaviour of the 3-point functions of scalar (" and 
tensor jij perturbations with correlation functions of the stress tensor of the dual QFT. By 
late times we mean the end of the holographic epoch, which should be the beginning of hot big 
bang cosmologylf] In subsequent sections we will use the hologragraphic formulae to obtain the 
cosmological predictions for a universe described holographically by a weakly coupled QFT. 

More precisely, we consider 3-point functions involving the curvature perturbation on 
uniform energy density slices and the transverse traceless tensor ■yij {■ya = and di%j = 0). 
These variables are defined such that in comoving gauge, where the inflaton perturbation 6ip 
vanishes, the spatial part of the perturbed metric reads 

gij = a^e^^\ei\ij = a^e^^(bij + ^ij + ^7ifc7fcj)- (3.1) 

(For fully gauge-invariant expressions to quadratic order see Section 2.2 of [B].) We find it 
useful to work in the helicity basis where 

= % (3.2) 

and the helicity tensors isi) are summarised in Appendix 1X1 

The holographic formulae for the in-in 3-point correlators are given by 

((C(gi)C(92)C(g3))) 



- 2(^((T(gi)T(g2,g3))) +cyclicperms.^ , (3.3) 



((C(gi)C('Z2)7^^^n93))) 
1 



^2 Jm[i?(g-i)]Im[i?(g2)]Im[A(g3) 



X Im 



((r(g-i)r(g-2)r(^«)(g-3))) - 2{Q^l''\m{'ix) + <dt\mm) 

-2((T(g-i,g-2)r(^3)(g-3))) -2((r(g-i)T(^3)(g-2,g-3))) -2((r(g-2)T(^^)(g-i,5-3)))l, (3.4) 



^In particular, we will assume a smooth transition to hot big bang cosmology. Developing a holographic 
theory of reheating is very interesting but will not be pursued here, see the comments at the end of Section 2 



of [5] for a preliminary discussion. 
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= -^(lm[B{qi)]lm[A{q2)]lm[A{q3)]'~' 



X Im 



((r(gi)r(^2)fe)r(^3)(9-3))) - ^(^(^2) + A{qs))e(''^^\q,) - i?(g-i)e(^^^^)(g-, 



- 2((T(g-i)T(«2^3)(^2,g3))) - 2((r(^^)fe)T(^3)(g-^,g-3))) _ 2{{T^^^\q,)T('^\qr,q2))) , 

(3.5) 

'llMAmy' X Im[2((r(«^)(gi)r(^2)(^2)T(^^)(g3))) - ^e(^i^^^«)(g,) ^(g,) 

i i 

- A(^{{T^-'^\q,)r^''-'^\q2,q3))) + CYdicpeims}j], (3.6) 

where, as noted previously, the double bracket notation indicates correlators with the momentum- 
conserving delta function removed, i.e., 

(C(gl)C(g2)C(^i)) = {27^fS{Y,maqi)aq2)C{q3))), (3.7) 

and similarly for other correlators. The coefficients A{qi) and B{qi) are related to 2-point 
function of the stress tensor (see (j2.3p ). while T and T^'*) are the trace and helicity-projected 
transverse traceless part ofTij, as defined in Appendix[Al Similarly, T, T^'') and T^*^''^) are the 
trace and helicity projections of the T tensor (j2.2p . The theta functions, \ 62*^ q{sis2)^ 
0{siS2) g^j-^jj @(siS2S3)^ represent specific contractions of the helicity tensors, and are described 
in Appendix El The imaginary part in these formulae is taken after making the analytic 
continuation 

N = -iN, qi = -iqi. (3.8) 

where N is the rank of the gauge group, see Section [2l and q = +\f^ is the magnitude of 
the momentum. 

The right-hand sides of (|3.3p - (|3.6p were obtained by first deriving, using standard gauge/gravity 
duality, the holographic 3-point functions of the stress tensor along general holographic RG 
flows, either asymptotically AdS or asymptotic to non-conformal brane backgroundqj and 
then analytically continuing to the cosmological case. These correlation functions are in a flat 



® The non-conformal brane backgrounds are asymptotically AdS in the dual frame [18]; the final formulae 
are however the same in both cases |4l|6]. 
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background and are defined as usual (for example) by the path integral formula 

(T,,,,(g-i)...r,„,„(g-„)) = / [ci^]T,,,,(g-i)...r,„,„(g„)e-^«^-[^l (3.9) 



where (p denotes collectively all fields of the boundary theory. We work with the Wick rotated 
QFT and correspondingly the bulk has Euclidean signature too. This is convenient because 
bulk regularity in the interior translates into the standard Bunch-Davies vacuum after analytic 
continuation to cosmology [Ij. The holographic relations (|3.3p - (|3.6p may also be viewed as 
computing the wavefunction of the universe extending the analysis of [19] (see also the recent 
[20 j ) to a general class of FRW spacetimes. 

Note that all terms appearing in the right-hand side numerators of (j3.3p - (j3.6p . except for 
the 3-point functions, are semi-local, i.e., in position space two of the operators are coincident 
(in momentum space this corresponds to the correlator being non-analytic in only one of the 
three momenta). Due to the non-analytic powers of momenta appearing in the denominators 
of these formulae, however, semi-local terms in the numerator generate contributions to the 
bispectra that are non-analytic in two momenta. Hence, as discussed in j4j in the case of the 
scalar bispectrum, they may contribute, for example, to 'local' type non-Gaussianity. Note 
that these terms may be computed unambiguously in perturbation theory (and we will do so in 
the next section). In contrast, ultra-local terms, i.e., terms where all operators are coincident 
(or equivalently terms analytic in all three momenta), are in general scheme dependent since 
their value can be changed by local finite counterterms. 

One often defines the correlation function of the stress tensor by first coupling the theory 
to a background metric gij^ differentiating w.r.t. gij and then setting gij = 6ij. One has to be 
careful, however, if one is to match with the unambiguous expression in p.9p . The 1-point 
function in the presence of the source gij is defined by 

m,{x)), = -^=-±- [ [d^]e-'<^^-^^-'^^^^ (3.10) 



Higher point functions are obtained by further functional differentiation and then setting 
gij = 5ij . This procedure leads to a new insertion of the stress tensor when functional derivative 
acts on Sqft['P'i dij], but also leads to additional semi- and ulta-local terms. One source of 



such terms are the factors of g{x) in (jS.lOp and for this reason some authors (see for 
example [13]) define the correlators without such factors, i.e., 



(3,11) 

9i3 —^ij 
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Note however that these correlators differ from (j3.9p (hence the different notation: < T... > 
instead of {T...)) because the stress tensor of the theory in a curved background also de- 
pends on gij and the functional differentiation leads to additional insertions of Tijki{xi,X2) = 
6Tij{xi) /6g^\x2)\Q- A careful evaluation of all such semi-local terms is given in Section 4.1 of 
[6]. This is the origin of most (but not all) semi-local terms in (|3.3p - ()3.6p . 

To connect with holography let us now recall that in gauge/gravity duality one identifies 
the fields parametrizing the boundary conditions for the bulk fields with the sources of the dual 
operators, and the (renormalised) on-shell gravitational action with the generating functional 
of correlation functions |2HI22j. In particular, the boundary metric 5(o)jj in Feffer man- Graham 
coordinates is identified with the source of stress tensor [23] . One could also envision a 
holographic mapping where the source of the stress tensor gij is related to the boundary 
metric g(o)ij via a local relation. Then, as was recently emphasized in [7], the two holographic 
maps would differ by (semi-) local terms. Given any such relation one can straightforwardly 
work out the corresponding holographic formulae, including all semi-local terms. One can 
(at least partially) fix this potential ambiguity in the holographic map by requiring that the 
(anomalous) Ward identities derived in gravity and in QFT match. For example, the matching 
of the Weyl anomaly in odd (bulk) dimensions |24t [25] fixes g^j = 5(o)ij • In even dimensions 
one would look to match the semi-local terms in the dilatation Ward identity of higher point 
functions. The formulae (I3.3p - (j3.6j) we use here were derived using the standard identification 
gij = g(o)ij and the definitions of ( and ^ij in (13. ip . 

4 Evaluating the holographic formulae 

In this section and throughout, we will work in Euclidean signature, and to leading order in 
gfgjj and 1/N. It will be useful to express our results in terms of the elementary symmetric 
polynomials of two and three variables, which we denote 

ai23 = qi + Q2 + qs, bi23 = qiq2 + q2q3 + q-m, C123 = 91^293, 

ai2 = qi + q2, bi2 = 9192, (4.1) 

with similar expressions for 023, 6231 stc. We also define 

= {qi + 92 + 93) (-91 + 92 + 93) (91 - 92 + 93) (91 + 92 - 93) 

= -ai23(a?23 - 401236123 + 8C123), (4.2) 
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and the quantities 

Pijkl = '2^i{k^l)j — ^ij^kl, Pijkl = ^i{k^l)j — Sij^kl- (4.3) 

Note that A is equal to 1/4 of the area of the triangle with length sides equal to the momenta 
QiiQ2,Q3 (Heron's formula). 

4.1 2-point functions 

In this subsection we recall the contribution to the 2-point function (j2.3p from each of the 
individual fields derived in [H [2J : 

A4, = B^ = ^M4,N^f, = ^M^N^q^ = 0, 

AA = BA = ^MAN'f, A^ = ^M^N'f, B^ = 0, (4.4) 



and thus in total we have 



where 

A/'{A) = AA4 + AA^ + AA^ + 2Af^, Af^B) =^A+ M^. (4.6) 

4.2 Contribution from minimal scalars 

The 3-point function for minimal scalars is given by the integral 

= N4,N PijabPklcdPmnef I [Ml ZoTT _ .n._ , _ (4.7) 

J q\q- qir{q + q2y 

Here, and throughout, we will make use of the shorthand notation [dq\ = d^g/(27r)^. Evaluating 
this integral using the general methods discussed in Appendix [Bl we obtain 

{{T^{qi)T^{q2)T'^\qz))) = - [Sauql + 2(3a?2 - 4612)^3 + ai2(3a?2 - 46i2)] , 

miqi)T^^\q2)T!p'\m = - o.^^tfr, {qi - -a23f[5ql + 20a23g1 + (29^23 + 6523)^1 

+ 023(170^3 + 21623)gt + 3ai3(«23 + 8^23)9? + 2a^3(a^3 + 3623)9? 
+ (3a23 - 6023623 - 326^3 )^i + 023(023 - 3623)] , 
11 



(ql - aia + 4^23)^ [5g? - (aia + 2623)91 + a23(ai3 - 3^2; 



8192 



2 \2 



32768V2 at,3c1,3 ^'''^^ " '"'^^'"^ + '"'^^'"^ " ' 



2 \2 



32768^/2 a223cf 23 



■ (93 - 012)^ [3g| + 4oi2gt + (a?2 " 2^12)9! 



+ ai2(a?2 - 36i2)gi + 4ai2(a?2 - 3612)93 + a?2(3af2 - 7612)] • 



(4.8) 

All remaining 3-point functions for minimal scalars may be found from these via permutations 
and/or a parity transformation. (The result for three insertions of the trace is given in |3].) 

Turning now to evaluate the semi-local terms in the holographic formulae, for minimal 
scalars 



Tf^.,,(xi, X2) = -\i5^,T^l + P^JklT^)S{S^ - X2). 

We thus have 

^Tt^{qi)'^tvrnniMz))) = -\hdTf^{qi)TL{-qi))) " \Pklmn{{Tf^{qi)TH-qi))) , 
from which we may extract the helicity-projected components 

{{mi)'^f\q2,m = -lB^{ql)Q^^'\q^l 
miqint^'\q2,q3))) = -B4q,)9(^^^-\q,), 

{{T^''\qi)T^{q2,m=0, 
((Tf^(9i)Tf^^ (92,9-3))) = -^A^(9-i)^(^^^^)(9-.), 
((rj^^^(9-i)T(^^^^^(9-2,9-3)))=0. 

4.3 Contribution from fermions 

In momentum space 



(4.9) 



(4.10) 



(4.11) 



9ym 



tr 



^Pijab I [dq]{qi - 2q)a : ^P^ {qjlb'ip^ {qi -q): 



(4.12) 



from which it follows that the 3-point function is given by the integral 
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-l\r Kr^p p p r f,^^-^ Qu{q-qi)v{q+q2)wi'2q-qi)ai'^q+q2)c{'^q-qi+q2)e 

— ^■'^tp 't^ijabJ'^klcd't'^rnnef^ ubvfwd / [ 



g2(g - gi)2(^ + 



(4.13) 



where 



^ubvfwd = iT^ilulblvlflwld) 

— '^^ub^vwdf ~l~ '^^uvPbwdf "^^ufPbwvd ~l~ '^^uwPbfvd '^^udPbfvwi (^-l^) 

recalling that {7i,7j} = —25ij. 

Evaluating the integral explicitly, we fincj^ 

((r^(g-i)r^fe)rj+)fe)))=o, 



((r^fe)r(+)(g-2)r(-)fe))) = -^-^ a23(a^3 - 3623)(g1 - + 4523)^ 



2 



2048 6; 



23 

2 \2 



{{Ti;\qi)T^-'\q2)T^-'\q3))) = - ,^„ftf_/_2 [a?23 " 2aL3"^i23 + a?23Ci23 + 32^3] , 

»iy2V^ai23Cl23 



n23H23 

r2\2 



((Ti+^(9-i)Ti+^(g-2)T(-)(,-3))) = - ,,,f"^I^ _2 - ai2)n^1 + ai2g1 - 612^1 

SiyzV^ ai23'^123 

+ a?2(a?2 - 3^12)^3 + a?2(a?2 - 2612)] • (4.15) 
The correlators with only one trace may be written in the condensed form 

m{qi)Ti:'\q2)Ti:'\q3))) = -l{A^{q2) + A4qs))e(^-'^\q.). (4.16) 

This result is in fact fully determined by the dilatation Ward identity (accounting for its 
semi- local nature), as we discuss in the next section. 



The result for the momentum dependence of the (+ + +) and (+ H ) components agrees with the 

resuh reported in 7,, but the overaU normahsation (after taking into account the difference in conventions, see 
Appendix |X| still differs from ours by a factor of four. Actually all correlators in 7 , including the 2-point 
functions, differ from ours by the same overall factor of four. 
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To compute the semi-local terms appearing in the holographic formulae, we find by explicit 
calculation that the operator 

T^-,,(^l,f2) = C^^^L^^nixi)Six, - X,) + C'^^,^J^{x,)dJ{x, - (4.17) 

where partial derivatives are taken with respect to xi, and the local operators 

Mmn = -^t^\^i^^lm%ni^^], Jm = ^tl[ji;^ J^tP^] , (4.18) 

9ym ^ 9ym ^ 

are associated with the coefficients 

^ijklm = ^i{k^l)j^mn + ^ij^m{k^l)n " ^ij^kl^mn - ^m{k^l){i^j)n- (4.19) 

As might be anticipated from their respective conformal dimensions (see Section [5]), 

{T^iM^n) = {TtiTD, {TtiJj = 0, (4.20) 

from which it follows that 

iiTfM)'^t^ni^2, qs))) = ci^L((^?(9i)7^at(-^-i)))- (4.21) 
Projecting into the helicity basis, the components appearing in the holographic formulae are 

m{qin'-;'''\q2,q3)))=0, 
{{Ti:'\.qi)r4'i2,m=0, 



{{Ti''\qi)T^;'\q2,m = -lA^{qM''''\q^), 
{{Ti:'\qi)Tl;'''\q2,m = -j^A^{ql)e^''''''\q^). (4.22) 
4.4 Contribution from conformal scalars 

As discussed in the introduction, the stress tensor xfj for minimal scalars may be decomposed 
as Tfj = xfj+Cij, where Tfj is the stress tensor for conformal scalars and Cij is an improvement 
term. For fields in the adjoint representation, and in momentum space, the improvement term 
takes the form 

Cijiqi) = fi^ijiqi) [ [dq] ■.(t>\q)<i>\qi - I) : 1 • (4.23) 
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Due to the presence of the projection operator vrjj, it follows that T^\q) = T'^\q) and hence 
the conformal scalar 3-point function involving three helicities is equal to that for minimal 
scalars. Similarly, the correlator 

{{mi)T^;'\q2)Ti^'\m = {{mi)Ti''\q2)Ti''\m - {mi)Tl''\q2)Ti''\m. (4.24) 

where the latter term may be evaluated from the integral 

mi)TtM2)TLm) = \M,mP,,.,P^^.^, [ [dg f"°(^\+^^)^(_^- ^^)^(g + g^)^ (4.25) 

^ J Q[Q-Qi){(l + (l2) 

Thus, to evaluate the conformal scalar 3-point function involving two helicities, only this 

integral needs to be computed since we already have the result for minimal scalars. Finally, 

evaluating the trace T^{q) directly, it is straightforward to show that the conformal scalar 

3-point function involving only one helicity vanishes. 

In light of these considerations, the 3-point functions for the conformal scalar field x arc 

((r^fe)r^fe)r(^^)fe))) = o, 

((r(^^)(g-i)T(«^)fe)r(^«)fe))) = ^ {{Tl''\q,)Tl''\q2)Tl'^\m- (4-26) 
Turning now to the semi-local terms in the holographic formulae, by direct calculation 



1 
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.(3) 



+ Ci;!,imnidmdn6ix,-X2))\OHx,), (4.27) 

where partial derivatives are again taken with respect to xi, the dimension one operator 

0>^ = ^tr[(x^)2], (4.28) 
9ym 

and the prefactors are 

*The result for momentum dependence of the (+ + +) and (+ H ) components agrees with the result 

reported in v2 of ;?], but the overall normalisation (after taking into account the difference in conventions, see 
Appendix |A]) still differs from ours by a factor of four, see also footnote [T] 
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^ijklmn ~ '^^ij^k(m.^n)l + h(k^l)j^mn — ^ij^kl^mn — '^Sm{iSj){k^l)n, 

^ijklmn = ^ijh{mK)l + Si{k^l)jSmn - ^ij^kl^mn - '^^m(i^ j){k^l)n + 5kl^i(m^n)j- (4.29) 

The precise form of these prefactors is not important, however, since 

{T^-O^) = (4.30) 
due to the differing conformal dimension of the two operators, hence 

{{nM^)nimn{M))) = -\5kdT^M^)T^n{-^l)))- (4-31) 

The hehcity-projected components appearing in the holographic formulae are then 

{{T(f^\q,)T^{q,,q,)))=0, 
{{Tif^\q,)T(^^\q2,q3))) = -^^^(g-i)0(^^^3)(^- )^ 
{{Ti^^\q,)T^^-^-^\q2,q3)))=0. (4.32) 
4.5 Contribution from gauge fields 

The stress tensor for the gauge fields is given by ()2.7p plus the contribution due to the ghosts 
and gauge-fixing terms. The latter contributions are BRST exact however, and thus they 
should not contribute to the 3-point function. Indeed, we show in Appendix [C] that their 
contribution cancels. 

Introducing the Hodge-dual field strength , where 

^ij ~ ^ijkGi, = — Pijab'^^ G^a^l, (4.33) 

9ym 

and evaluating its propagator, one finds the contribution from gauge fields to the 3-point 
function is given by 

((7;j(gi)rfc4(g2)r^„(g3))) = -J^AN^PijabPklcdPmnej j WWc{q)T^be{q ' qi)^^ df{q + q2) ■ (4.34) 

Upon closer examination, this integral may equivalently be expressed in terms of the 2- and 
3-point functions for minimal scalars. 
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^((7;j(g-l)T,4fe)Ti„(g-3))) = ((i;t(g-l)T,*,(g-2)T^„(g-3)>> -g./™na6((T;t(g-l)r,t(-g-l))) 

- Qmnijab{{Tti{q2)T^b{-q2))) - Qzjklab{{Tin{Q3)T^bi-q3))) , 

(4.35) 

where 

Qijklmn — PijacPklbcPabmn- (4.36) 

This result is a consequence of the fact that Gj may be identified with the operator 9j(/>, where 
(p is a massless scalar field. The appearance of the various semi-local terms in (j4.35p then 
reflects the fact that diGi vanishes identically, while d'^cj) vanishes on-shell only. 

Turning now to evaluate the semi-local terms appearing in the holographic formulae, a 
short calculation reveals the operator 

1 



'^ijkl{^1^^2) = -^^SijT^i + PijklT^ + QijklmnT^n\H^l " ^2)- (4.37) 

Making use of the fact that the 2-point functions for gauge fields and for minimal scalars 
coincide (see (|4.4|) 1. it then follows that 

^((i;f(g-i)T^,_(g-2,g-3))> = {{Tt^{qi)Ti^Jq2,m - \Qklmnab{{Tf^{qi)Ti{-qi))) . (4.38) 
Thus, from (j4.35p . 

^ [iT,j(q,)T^l{q2)T^^m) " 2 (((T^f (gOT^.^Jg^, 93))) + cyclic permsj 

= {{Tf^{qi)TtM2)TLm) - 2(((T;t(^-i)TLnfe,9-3))) +cyclic perms). (4.39) 

This particular combination, suitably projected, appears in the numerator of all the holo- 
graphic formulae for cosmological 3-point functions (namely (13. 4p . ()3.5p and (13. 6p ). Thus, 
since the 2-point functions for gauge fields and minimal scalars also coincide, we see that 
gauge fields and minimal scalars necessarily make identical contributions to all cosmological 
3-point functions. Since scalars and vectors are dual in three dimensions, this result is per- 
haps not unexpected, and indeed similar behaviour was noted in [1] for the case of the scalar 
bispectrum. 



iA I I') rriA I /■ \ rjiA 
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5 Ward identities 



In the previous section we computed all relevant 3-point functions and semi-local terms by 
direct computation of 1-loop Feynman diagrams. In this section we elucidate the structure of 
these correlators by ascertaining the extent to which they are determined by Ward identities. 

5.1 Minimal scalar s from conformal scalar s 

As noted previously, the stress tensor for minimal scalars may be decomposed as 



Ox 



2 -^Ai^ 

Sym 



J\2i 



(5.1) 



where T/- is the stress tensor for a conformal scalar field and Oi is a dimension one scalar 



operator. The 3-point functions of Xf- may thus be expressed in terms of 3-point functions of 



the conformal fields T^j and Oi. Specifically, we find 



.(^3) 



-.(ss). 



=;(S3)/ 



(S3), 



=i(s3)/ 



+ |((r^(g-i)Oife)r(^^)(g-3)» + ^|((Oi(gi)Oi(g-2)r(^^)fe))), 



{{T^{q,)T^{q2)T^{qs))) = ((T^(^i)r^fe)r^fe))) + 



((Oi(gi)r^(g2)r<^fe))) + 2 perm. 



+ 
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((01(51)01(^2)^4^3))) + 2 perm. 



+ ^((Oi(g-i)Oi(g-2)Oi(g-3))). 



(5.2) 



Recalling that gauge fields contribute the same as minimal scalars, the computation of general 
3-point functions thus reduces to computing a set of 2- and 3-point functions in a CFT. (Note 
that free fermions are also a CFT.) These correlators are in turn (almost) uniquely determined 
by Ward identities, as we now show. 



5.2 Trace Ward identity 

In light of the above, we are interested in correlation functions of the stress tensor and of a 
scalar operator Oa of dimension A = 1 in a three-dimensional CFT. We will be more general 
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however and discuss the case of any d and A, provided only that A ^ d. The trace Ward 
identity in the presence of a source (po for O/^ reads 

{T{x))s = {A-d)MOA{x))s. (5.3) 

This Ward identity imphes that n-point functions involving an insertion of the trace of the 
stress tensor are given by semi- local terms involving (n — l)-point functions. These relations 
can be obtained by functionally differentiating (j5.3p w.r.t. the sources (n — 1) times and then 
setting them to zero. Noting that all 1-point functions vanish, for the 2-point functions we 
find 

{T{x)OA{y)) = 0, {T{x)Tki{y)) = 0, (5.4) 
and for the 3-point functions, 

(T(x)Oa(2/)Oa(^)) = {§^,Oa{z)) + (-g^OA(y)) 

+ (d-A) [dix - y){OA{x)OA{z)) + 6{x - z)(Oa(x)Oa(2/))] , (5.5) 

{nx)TMOAiz)) = 2{^^Oa{z)) + i^TM) + (T(x)g^) 

og'^^iy) o<po{z) 000 (^) 

+ {d-A)6{x-z){ni{y)OA{x)), (5.6) 

{nx)TMT^M)) = 2(^r_(.)) + 2(^p^r,,(y)) + 2(r(x)pM). (5.7) 

As the stress tensor in the presence of sources has (po dependence 

Tij [g, 0o] = Tij [g, 00 = 0] - fify^oCA, (5.8) 

we may in addition identify 

^^ = -d6ix-y)OAix). (5.9) 

Then, since the CFT correlator {TijOA) vanishes for any operator Oa with dimension different 
to d, equations ()5.5p and ()5.6p reduce to 

{T{x)OA{y)OA{z)) = -A[5{x-y){OA{x)OA{z)) + 5{x-z){OA{x)OA{y))\, (5.10) 

{T{x)ni{y)OA{z)) = 2{J^^OAiz)). (5.11) 

Finally, it is convenient to express (j5.7p . (j5.10p and (j5.1ip in momentum space in terms of 
the T tensor defined in ()2.2p . Projecting into the helicity basis, we obtain the complete set of 
trace Ward identities 

{{T{qi)OA{q2)OAm) = -A[{{OA{q2)OA{-q2))) + ((OA(g3)OA(-^3)))] , 
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{{T{q,)T{q2)OA{m = 2((T(gi, g2)OA(g3))), 
{{T{qi)T(^^\q2)OAm) = 2{{T^'^\qi, q2)OAm) , 

{{T{q^)T{q2)Tm) = 2[{{T{q,)T{q2, qs))) + {{T{q2)r{q,,q^))) + {{T {qs)T {q^ , q2)))] , 
((r(g-i)T(g-2)T(«3)(^-3))) = 2[{{T{q^)T(■':^\q2,q3))) + {{Tiq2)T^'^^Hqi, qs))) 

+ {{riqi,q2)T('^-^\q3)))], 
{{T{qi)T^^^\q2)T('-^\m = l{MQ2)+Aiq3))^^^^^^ 

+ ((r(^^)fe)T(^^)(g-i,g-3))) + {{T^''Hq3)T^''\qi,q2)))], (5.12) 

where A{q) is the transverse traceless piece of the stress tensor 2-point function defined in 
(j2.3p (for conformal fields the trace piece B{q) vanishes as a consequence of (|5.4p ). 

Comparing with our holographic formulae (j3.3p . (|3.4p and (|3.5p . we immediately see that 
conformal fields make no contribution to the numerators of these formulae, as found earlier 
by explicit calculation. An important consequence of this, as we will see in Section [6l is that 
the CCC) CC7 and ("77 cosmological shape functions are forced to be independent of the field 
content of the dual QFT. 

Further insight may be distilled from the trace Ward identities (I5.12p by replacing the 
semi-local contact terms on the r.h.s. with 2-point functions of Tij and Oa. On general 
grounds, the T tensor has an expansion in terms of local operators of dimension less than 
or equal to d, and for fermions and conformal scalars we computed this explicitly in ()4.17p 
and (j4.27p . Then, as we found in the analysis leading to (|4.2ip and (|4.3ip . only operators of 
dimension d contribute to the correlator {TijTkimn)^ permitting it to be expressed in terms 
of (TijTi^i). Substituting into (|5.12p our previous results (|4.22p and (j4.32p for the semi-local 
terms {TijTkimn), we obtain 

{{T^{q,)T^{q2)T(f^Hm = {{T4'li)TA'l2)T;;'\m = 0, 

((r^(g-i)r(^^)(^-2)r(^^)(g-3))) = -\{A^{q2) + A^{qs))9^^-^'Hq^), 

{{T^{qi)T^''\q2)T^''\m = -\{A^{q2) + A4qs))e^^-^-^\qi). (5.13) 

Thus, all our earlier results in (|4.15p involving the trace are in fact a consequence of the 
Ward identities (noting also ()4.16p ). and similarly for all our results in (I4.26P involving T^. 
For the latter, note that A^{q) = {M^/2M^)A,^{q) from (j4.4p . hence from ()5.13p we have 

((T^(g-i)T(^^)(g-2)r(^«)(g-3))) = ^((T^(ft)Ti^^^(g2)ri^^^(g-3)». (5.14) 
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As well as confirming earlier calculations, these formulae additionally serve as a check of the 
overall sign in our 3-point function integrals. 

To check the results of our 3-point function calculations for minimal scalars using (j5.2p . 
we must also evaluate the semi-local terms on the r.h.s. of (j5.12p involving the correlator 
(CiT^.^^), where T^.^^ denotes the T tensor for conformal scalars. The expansion for this 
latter quantity may be read off from (|4.27p (replacing x with 0). The correlator (OiT^.^^) 
receives contributions only from terms of dimension one in this expansion, and so we find 

{{^I;'\qi,q2)0,{qs))) = l-f^Qi^-\q,){{0,iqs)0,{-qs))). (5.15) 
Substituting these expressions into (|5.12p . we obtain 

{{f4q^)Oi{q2)Oim) = -((Oi(g2)Oi(-g2)» - ((Oi(g3)Oi(-93)», 
((r^(g-i)r^(g-2)Oi(g-3))) = ^(a?2 -26i2 -g1)((Oi(g-3)Oi(-g-3))), 
{{f4q^)f^''\q2)0,m) = lg-|e(^^)(g-.)((Oi(g-3)Oi(-g-3))), (5.16) 
{{T^{q{)T^{q2)f^{m = 0, (5.17) 

where for the latter equation we used (I4.3ip . The trace Ward identities thus supply all terms 
appearing on the r.h.s. of ()5.2p that involve the trace T^. 

5.3 Conformal Ward identities 

In the previous subsection we showed how the trace Ward identities determine the 3-point 
functions involving the trace of the stress tensor in terms of 2-point functions. Thus to 
determine all correlation functions, it remains to obtain 

(T^s,)j,{s2)r^{s,)^^ (OiTJ'^^tJ"^)), (OiOiTJ'^^), (OiOiOi), (OiOi). (5.18) 

These may be directly computed using the methods described in Appendix [Bj For confor- 
mal scalars and fermions, the result for (T^^^^T^^^^T^^^^) is given in (j4.15p and (j4.26p . The 
remaining correlators are found to be (we suppress a common overall factor of iV^A/^ ) 

{{0,{q)0,{-m = 

iq 

((Oi(gi)Oi(g2)Oi(9-3))) = 3^, 

Cl23 

21 



(2gi + a23) 
I6V2 al^^hsqi 



((rj+)(g-i)Tj-)(g-2)Oi(^-3)» = ^^^^^ (^i - + 46i2)'(a?2 + 2bu - 5^1), 
rr^fe)Ti^^(^-2)Oi(9-3)» = ti'- [ - - 20^1291 - (29a?2 + 6612)^1 

2U48 ai230i293 

- 8012(20^2 + 36i2)g| + a?2(^i2 ~ 36612)^1 

+ 4a?2(a?2 - 6612)^3 + a% " 60^2^12 + 326?2] • (5-19) 



Subsituting these expressions into (j5.2p . along with those in (|5.16p . we recover ah the results 
for minimal scalars listed in ()4.8p that involve the trace T^. (One may additionally check we 
recover the result for (T^T^T^) in equation (101) of [4J.) 

The correlation functions (|5.19p are almost uniquely determined by conformal Ward iden- 
tities. This was discussed in position space in |14j . and more recently in momentum space 
in [7] (see also J26j). More precisely, the 3-point function of the stress tensor is unique, up 
to two constants (which is our case may be taken to be the number of conformal scalars and 
the number of free fermions), and all remaining 3-point functions are unique up to an overall 
constant. Thus, a non-trivial check of the correlators listed above is to verify that they indeed 
satisfy the special conformal Ward identities. (The scale Ward identity is satisfied by inspec- 
tion.) These identities take the form of differential equations that the correlators must satisfy, 
and are listed explicitly in Appendix lD.il (they may be obtained by Fourier transforming the 
position space Ward identities whose derivation is discussed, for example, in [27]). We have 
checked that our conformal correlators satisfy these Ward identities. 

6 Holographic predictions for cosmological 3-point functions 

Having computed all relevant QFT quantities we can now evaluate the holographic formulae. 
It is instructive to first use the trace Ward identities and the relation of minimal scalars 
to conformal scalars in order to express the cosmological 3-point functions in terms of CFT 
correlations functions. This yields 





X 



[qlqlql{{Oi{qi)OMOM))) 

+ {2ql{ql - ql - ql){{Oi{qi)Oi{-qi))) + 2 perm.)] 
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iqU0i{qi)0^{q2)f^''\qs))) 
+ [qfe[''\qMOiiqi)Oii-qi))) + (q, ^ q^) 



,14 



N^Ml^^qlqUAO,{q,)0,{-q,))) 

\Qql{{0,{qr)fl''\q2)fl'^\q,))) 
+ q\ f2e(^2^3)(^^)_e(.2.3)(^.) 



{{aqi)i'^''\q2)i^''\qz))) 



((Oi(<7i)Oi(-gi)))J , (6.1) 

where ■N'^a) -^{b) ^I'e defined in (14. 6p and we have made the dependence on the number 
of fields exphcit by considering the Oi and T^'^^ correlators to be those of a single field. (The 
analytic continuation p.Sp has also been implicitly performed; the correlators appearing above 
are therefore those in ()5.19p with qi replaced by qi ) As discussed in the previous section, the 
trace Ward identities imply that the numerators of the holographic formulae (j3.3p - (|3.5p for the 
above correlators receive no contribution from conformal fields and are therefore proportional 
to A/'(B), the number of non-conformal fields. The dependence of these correlators on the 
field content is then simply given by an overall factor, amounting to ■N'{b) divided by the 
corresponding factors in the denominators of the holographic formulae. 
For the 777 correlator, we find 



((^(^i)(gi)^(^2)(^2)^(s3)(^3))) 



,24 



N^^fA)^lqlql 



(S3)\ 




(6.2) 



considering again the correlators to be those of a single field so as to make the dependence on 
the number of fields explicit. It then turns out that 

e(++-)(g.) ' 



1024 



(6.3) 



1=1 



i.e., these correlators differ only by the helicity projection of a semi-local term. Thus, while the 
3-point function of the stress tensor at separated points in general depends on two constants 
in d = 3 |14] . only one combination survives the (+ H — ) helicity projection. (This was also 
shown in [7] using the conformal Ward identities.) The specific form of the semi-local term 
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in (j6.3p is then such that the ^(+)^(+)^( ) correlation function depends on the field content 
through an overall multiplicative constant only. On the other hand, 

A2 



(+)\ 



//7.(+)'r(+)'r(+)\ 



+ — — — 



1024 



+ 



Cl23 



* 128^/2 a'' 



(6.4) 



123 



and so these correlators differ by the helicity projections of both a semi-local and a non-local 
term. The non-local term reflects the fact that both solutions for the 3-point function of the 
stress tensor at separated points survive the (-|- -|- -|-) helicity projection, and leads in turn to 
a more complicated dependence on the QFT field content in the ^(+)^(+)^(+) correlator. 

Returning to (j6.ip and substituting in our results for the remaining correlators, we first 
recover the result derived in [1], 



((C('Zi)Cfe)Cfe))) 



512 



512 



i 

A2 



-3 



-2qiq2q3 - ^Qi + {(I1Q2 + 5 perms)) 



^flB)N^ ai23cf23^ 



(6.5) 



showing that the scalar bispectrum exactly coincides with the equilateral template. In the first 
line, note that all the terms but the one proportional to qiq2Q3 originate from semi-local terms 
in the numerator of the holographic formula (j3.3p . Without their contribution we would not 
have been able to distinguish the equilateral shape from others involving a similar factor of 
Q1Q2Q3 in the numerator (for example, the orthogonal shape [28], for which the corresponding 
numerator is —8515293 — ^^^gf -|- 3{qiq2 + 5 perms)). In fact, due to the factor of Hi ^i"^ 
coming from the product of 2-point functions in the denominator of the holographic formula, 
the semi-local term Yli Qi in the numerator generates a contribution to the bispectrum of 
exactly the 'local' type. It is therefore essential to include the contribution of all semi-local 
terms in the holographic formulae, as we have been careful to do. 
For the remaining correlators, we find 

2048 A2 



((C(9l)C('Z2)7^+^('Z3))) 

((C(9i)7^+H'?2)7(+)(93))) 

((^(+)(^l)^{+)(g2)7(+)(53))) 



'^123'^12393 



(^123 ~ 0^123^123 — C123) — 012353 



512 



(Qi - alsf ili - 43 + 2^23) + 



32bi 



23 



'-123 



512 



1024 



^23 



23A 



A ai23 



'123 



('^123 



ai23^ 



123 



C123J 



(A) 



64c|23- 
"123 ■ 
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1024 

{{¥'^Hqi)l^^\q2)l'^'\q3))) = /K.r4^^2 2 —(is - «12)''(a?23 - ai23&123 - C123). (6.6) 

V ^iV A/^^^ "123<'123 

We would now like to define corresponding shape functions, i.e., we wish to write these 
correlators as bispectra: a product of power spectra times a shape function. To do so, we first 
define the dimensionless 2-point amplitudes 

^(CC) = qHiCiM-q))) = ^^2^, ^(77) = q'{{i^^\qn^^\-q))) = (6-7) 

and similarly the dimensionless 3-point amplitudes, e.g., 

Am^+^) = ghhs ((C(9i)C(g2)7(+H93))), (6.8) 

with analogous expressions for the other correlators. Physically, these quantitites parametrise 
the contribution per logarithmic interval of wavenumbers to the corresponding position-space 
expectation values with all insertions at the same point, e.g., 

^^'^^^^ = (C'(i^)¥+H^)) = ^/ (lldlnq}j ^(CC¥+^), (6.9) 

where the latter integral ranges over all possible triangle side lengths in momentum space. (For 
reference, the usual logarithmic power spectrum is simply = (l/27r^)„4(CC)-) The dimen- 
sionless 3-point amplitudes may now be naturally re-expressed as a product of dimensionless 
2-point amplitudes and a purely momentum-dependent shape function: 

^(CC7^'^^) = ^(CC)-4(77)5(CC7^'^^ ^(C^'^H^^^^) = ^'(77)5(C7^'^H^'^^), 

^(^(5i)^(s2)^W) ^ A'^{J^)s{^^'^^^'^^^'^). (6.10) 

Explicitly, the shape functions are given by 

1 r 

^(((7^"^^) = TlK^ 2 ("123 - «123&123 - C123) - 0123^3 



326i3 



SiC^MjM) = - ahf [(,? - ah + 2623) + 

S{Cj(+)j(-)) = --^{q! - al, + 4623)^(9? - aia + 2623), 



4 

123 



128613 



5(^W^(+)^W) = [(a?,3 - ai236i23 - C123) - (l - 4 

64v2 C123 ^ 



A/",/; \ 64c' 



3 



(A)' ^123 
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= 1 _^(g3 _ ai2)"(a?23 - ai236l23 - C123). (6.11) 
b4:^/2 0123^123 



Thus, with the sole exception of 5(7('''^7(+^7*^^)), all the shape functions defined in this manner 
are independent of the field content of the dual QFT. (Indeed, this was our motivation in 
selecting the factors of A{CC) and ^(77) appearing in (j6.10p .) Prom our previous discussion, 
we see that for the shape functions involving one or more factors of C this property is a 
consequence of the trace Ward identities, which limit the field content-dependence of the 
corresponding bispectra to a single overall factor. The independence of 5(7*^+^7(+)7^^-') from 
the QFT field content arises similarly from the fact that the corresponding bispectrum depends 
on the field content via an overall factor only. As we saw above, this latter property relies on 
both the conformal Ward identities and the precise form of the semi-local terms appearing in 
the holographic formula. 

We have plotted the holographic shape functions in Figs. [2] and [3] (along with their coun- 
terparts for slow-roll inflation which we discuss in the next section). In these figures we have 
adopted the expedient of scaling all momenta such that qi + q2 + Q3 = ^ (note that the shape 
functions are invariant under a constant rescaling of all momenta). By the usual triangle 
inequalities, the allowed range for any two momenta, say qi and q2, is then < qi < 1/2 and 
1/2 — gi < q2 < 1/2 as displayed. In each case, we have chosen to plot the two momenta under 
whose interchange the shape function is symmetric. 

Note that the usual plotting convention adopted for the scalar bispectrum S{(CC) (namely, 
ordering the momenta > 52 ^ qs and then scaling qi to unity, with the triangle inequality 
then constraining q2 > 1 — 93) is not applicable to the correlators considered here, since in each 
case (with the sole exception of 5(7^^)7('''^7(^^)) one of the three momenta is distinguished 
and so the required ordering of momenta cannot be accomplished without loss of generality. 
Without this initial ordering step, rescaling one of the momenta to unity then fails to yield an 
upper bound on the magnitude of the remaining momenta, resulting in a plot with unbounded 
area. This problem is neatly sidestepped by constraining the total perimeter of the triangle 
to be unity, instead of the length of one the sides. 

7 Comparison with slow-roll results 

Slow-roll inflation predicts the correlators of three gravitons are 

64\/2 ^123 



m^^Hqi)l^^Hq2n^^HQ3)))sR = ^ ^^^(«?23 - ai23fci23 - C123) 
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(a) 5(7(+>7(+'7*"^) 




(e) 5(CC7'+') (f) 5s«(CC7<+') (g) A5(CC7'+') 



Figure 2: Isoperimetric plots displaying the holographic and slow-roll shape functions, 
as well as the difference between them {e.g., A5(7^'^^7^^^7'-^-*) = 5(7'''*''*7'''^-*7''"'''*) — 

invariance of the shape functions under a rescaling qi — t- Xqi of 
all momenta has been exploited to set qi + q2 + Qs = li constraining the allowed momen- 
tum values to those displayed. Each plot is symmetric under interchange of the appropriate 
momenta as expected. Note that 5(7("^)7(~'^^7(~-') (shown in plot (a)) coincides for the holo- 
graphic and slow-roll models. In plots ()2bp and ()2dp we have set = J^(A) to maximise 
A5(7(+)7*^+)7*^+)) for illustrative purposes. 
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(a) 5(C7<+'7*+') (b) 5sii(C7'+H^+') (c) A5(C7<+H'+') 




(d) 5(C7*+H*"') (e) 5sfl(C7'+'7'"^) (f) A5(C7'+'7*"^) 



Figure 3: Isoperimetric plots for holographic and slow-roll shape functions continued. In plots 
()3dp and (i3ell note that both shape functions are actually finite along the line q^, = 1/2 — ^2 
(i.e., = 1/2); we have simply restricted the plot range to exhibit the overall shape more 
clearly. 
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{{¥^Hqi)l^^Hq2W'\q3)))sR = -—7= — 5— (93 - «12)^(a?23 - 0123^123 " Cm)- (7.1) 

b4:^/2 0123^123 



We may recover these results exactly from our holographic model by setting 

2AA^ = Af^+Af^+M^, ^N'Af^A) = J^. (7.2) 



In particular, the latter relation is also consistent with matching the amplitude of the graviton 
2-point function of slow-roll inflation and the holographic model. The first relation is consistent 
with that found in v2 of [7] for the special case where Ma = Mcf, = 0. (Note however that our 
careful treatment of the semi-local terms in the holographic formulae enables us to correctly 
recover the entire slow-roll bispectrum (j7.ip .) For general QFT field content (for which the 
first relation in ()7.2p is not satisfied), the and ■y^^^ ^^^^ ^^^^ holographic bispectra 

in (|6.5|) coincide precisely with the corresponding bispectra derived in [7j for slow-roll inflation 
in which one includes an additional term in the action proportional to the Weyl tensor cubed. 
Relative to [7J, our QFT additionally contains non-conformal fields, and our treatment of the 
semi-local terms enables us to recover the cosmological result exactly. 
The remaining slow-roll results are 

{{C{qi)C{q2)¥^Hq3)))sR = -—J=- 2 3 2 hl23 " ai236l23 - C123] , 

ibv^e* a'123'^123^3 

l^'^H^ , r) \o\ / 9 . „, \ ^^23 



{{aqi)i^'-\q2)i^'-\qz)))sR = -1^55^(9? - "23)' [{ql - 4, + 2623 



gi«123 



1^4 174 o/,2 

rv 1± ^ / 9 O \0r/9 9 ^-1 \ Ol/q 

128 6|3gf ' - ■ ^— - 910123' 

(7.3) 



{{aqi)l^^Hq2)l'--Hq3)))sR = -T^^TT^iq! - als + 4623)' M " a',, + 2623^ 



While these differ from the predictions of the holographic model, interestingly the difference 
is only in the last term. 

Evaluating the shape functions, for slow-roll inflation the 2-point amplitudes defined anal- 
ogously to ()6.7p are 

Asr{CO = ^, AsRm) = ^'Hl (7.4) 

The slow-roll shape functions then differ from their holographic counterparts by at most a 
single term: 



4V^ai23Cl23 ' 
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4ci23 

3 

123 



16ai23Ci23 ^ a' 

16ai23Ci23 

5sk(7^+H^+H^"^) =5(7^7^+^^"^)- (7.5) 

The holographic and slow-roll shape functions, as well as the difference terms in the expressions 
above, are plotted in Figs. [2] and [3l From these figures it is apparent that the holographic 
and slow-roll shape functions share the same broad quantitative features in all cases except 
for here, Ssr{C^^^^'J^^'^) has a simple pole as the momentum qi associated with ^ 

vanishes, whereas the corresponding holographic shape function has a zero. 

At a more quantitative level, a rough indication of the distinguishability of the holographic 
and slow-roll shape functions may be obtained by evaluating the cosine orthogonality measure 
proposed in |29j (following earlier work in [30]), 

C(S,S') = ^^^d^=, (7.6) 

where the weighted inner product 

F{S,S') = / dgidg2dg3^— '5(gi,g2,g3)'5'(gi,g2,93)- (7.7) 
J ai23 

Writing q\ = agi, q2 = aq2 and ^3 = a(l — qi — (72), the integral over a in the inner product 
factors out, since all shape functions we consider here are scale- invariant, i.e., independent 
of Q. This overall factor may then be discarded since its contribution to the cosine measure 
C{S,S') cancels between numerator and denominator. We may thus replace (|7.7p with the 
two-dimensional integral 



F{S,S') = j d^idg25(gi,g2, 1 - gi - q2)S'{qi,q2, 1 - gi - 92), 



(7.8) 



where the shape functions here are precisely the isoperimetric shape functions plotted in Figs. [2] 
and[3l 

Naively, one might expect the domain of integration would be < gi < 1/2 and 1/2 — gi < 
q2 < 1/2. Since however several of the shape functions have poles when one or more of 
the triangle sides are taken to zero, as we see from Figs. [2] and [3l one must further restrict 
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the domain of integration in order to obtain finite inner products. The physical justification 
for tliis procedure is that any real observation is only sensitive to momenta in some range 
Qmin < Qi < Qma.x- We will therefore restrict all rescaled momenta cji > e, where the cutoff 
e = 9min/2(Zmax ~ 5 X 10^'*. The domain of integration < gi < 1/2 and 1/2 — qi<q2<l/2 
is thus further restricted by the conditions qi > e, q2 > e, and I — q2 — Qs > ^- For shape 
functions with poles at the corners, the orthogonality measure (j7.6p will depend on the cutoff 
e, reflecting the fact that our ability to resolve the shape functions concerned depends on how 
sensitive we are to the corners of the distribution. 

Having thus carefully defined the orthogonality measure, one may now numerically eval- 
uate the orthogonality measure between each holographic shape function and its slow-roll 
counterpart. Rounding to two decimal places, 

^(^(+)^(+)^(+)) ^ 1^00, C{(:^^+^^+'>) = 0.33, C7(C7^+H^"^) = 0.67, C{CC¥^^) ~ 1-00. 

Values close to unity indicate nearly indistinguishable shape functions, while smaller values 
correspond to shape functions that are more orthogonal. (For comparison, the overlap between 
the standard local and equilateral shape functions evaluates to C = 0.34 with our cutoff 
prescription.) Overall, these values confirm one's impression by eye from Figs. [2] and [3l 
namely, that the holographic and slow-roll shape functions are nearly indistiguishable for the 
cases ^(+)^(+)^(+) and CCl^^\ while in the case ^^(+)^(+) the two shape functions may be 
distinguished by the presence or absence of a pole as the momentum qi associated with 
vanishes. 



8 Discussion 

In this paper we computed the complete set of bispectra (and defined and extracted the 
corresponding shapeqj) for a class of holographic models of the very early universe based on 
perturbative QFT. The leading 1-loop result actually depends only on the free part of the 
QFT, so in particular our results are also the complete answer when the dual QFT is free. 
The field content of the dual theory includes gauge fields, massless fermions, massless minimal 
and conformal scalars and thus the parameters that can appear in the results are the number 
of species for each type of field. The bispectra could, a priori, depend on these in a complicated 



® To our knowledge, shapes other than those relevant for purely scalar or purely tensor bispectra have not 
been discussed before. 
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way, but it turns out that we get instead (nearly) universal results that are independent of 
all details of the dual QFT, within the class of the theories we consider. Thus, these models 
make clean and precise predictions. 

One can trace this universality to the specific form of the holographic map, the fact that to 
leading order the QFT is free, symmetry considerations and properties of d = 3 theories. Let 
us explain this. Firstly, in three dimensions, vectors are dual to scalars so one may anticipate 
that the contribution due to gauge fields (at 1-loop order) is equal to that of the contribution 
due to mininal scalars, and we indeed find this to be the case. Taking this into account, the 
answer could then depend on three parameters, the number of conformal scalars, 7V^, the 
number of fermions, A/"^ and the total number of gauge fields plus minimal scalars, M(^b)- The 
trace Ward identity of the dual QFT and the specific form of the holographic formulae then 
imply that, in all correlators involving at least one factor of the field content appears only 
as a multiplicative factor and is such that the corresponding shape functions are completely 
independent of the field content. 

Let us now turn to correlators involving only tensors: these are effectively determined by 
the 3-point function of the stress tensor of a CFT. In three dimensions, this 3-point function is 
parametrised by two constants, which in our case are related to the field content. Indeed, the 
shape corresponding to three positive helicity gravitons docs depend on the field content, but 
surprisingly the shape for two positive and one negative helicity graviton is independent of the 
field content. This can be explained in part by the fact that the (T^+^T^+^T^"^) correlator (at 
separated points) is actually uniquely fixed by conformal invariance up to a single constant. 
We emphasize however that this by itself is not sufficient to explain the independence of the 
corresponding shape function from the field content, as the specific form of the semi-local 
terms (both in the holographic map and in (T^+^T^+^T^"))) is crucial for this to happen. 

Our calculations carefully include all such semi-local contributions. In the holographic 
formulae for the bispectra, these contributions appear as terms in the numerator that are 
non-analytic in only one of the three momenta. Since the denominator of the holographic 
formulae is however non-analytic in all three momenta, the net contribution of these semi-local 
terms to the bispectra is in fact non-analytic in two of the three momenta. Semi-local terms 
in the holographic formulae may thus contribute, for example, to 'local'-type non-Gaussianity 
behaving as l/(?i(?2 + perms. Contributions of this nature therefore play a crucial role in 
allowing different cosmological shapes to be distinguished. 

To get a feeling for our results we also computed the corresponding slow-roll results and 



32 



compared them with the holographic results. Firstly, comparing the power spectra one obtains 
a relation between the parameters N'^ , J^(a) ^{b) of the QFT and the parameters , 
and e^, of the slow-roll model. Comparing the 3-point functions, we find that the ^(+)^(+)^(~) 
correlators agree exactly, while the ^(+)^(+)^(+) correlators can be made to agree if one imposes 
that that the field content satisfies the relation 2A/',/, = M(f,+MA+J^x- explained in [7], these 
slow-roll correlators are constrained by the late-time de Sitter isometrics to satisfy conformal 
Ward identities, and thus at separated points they should be expressible in terms of the 3- 
point functions of conformal scalars and free fermions. Indeed, the linear combination found 
in v2 of [7j is the same as the one we find (setting A/"^ = Ma = in our relation). By taking 
into account the contribution from semi-local terms, however, we are further able to correctly 
recover every individual term appearing in the graviton bispectra. 

There is no apparent reason for the remaining slow-roll and holographic correlators to 
agree. Nevertheless we find rather similar results. To quantify the difference we used the 
cosine orthogonality measure of [29] to obtain a first indication of the distinguishability of the 
corresponding shapes. We find that the shapes for CCt'''*'' ^"^^ nearly indistinguishable, while 
for ^^(+)^(+), the two shapes may be distinguished (as a consequence of differing behaviour 
in the squeezed limit where the momentum associated with the Q goes to zero), with the case 
of lying in between. 

All in all, we have a rather complete understanding of this class of models and their 
phenomenology. There are still a few things to be understood better: what constrains the 
semi- local contributions to the tensor correlators, and why are the holographic results appar- 
ently close to slow-roll ones? One can presumably also understand the squeezed limit of the 
correlators using Ward identities. On a whole, however, the structure of these models is rea- 
sonably firmly understood. It would be interesting to arrive at a similar level of understanding 
for the class holographic models that are based on deformations of conformal field theories. 
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A Helicity tensors 

This appendix summarises our notation and conventions for heUcity tensors and their con- 
tractions. To facihtate the comparison of our results with those of [7], we also briefly review 
the spinor helicity formalism of this latter work. 

We use helicity tensors e -^^ (q) satisfying the standard identities 

^^M^) = leffiHfi-^, = 2^^^'- (A.l) 

where helicities Si take values ±1 and our conventions for are those of f3T] (see p. 233)). We 
may go from a tensor basis to a helicity basis by contracting with e^^"* (q) . Explicitly, the trace 
and helicity components of the stress tensor Tij and the Tj^fc; tensor are defined by 

r(^) = 6,,T,,iq), rW(g-) = ie(;)(-^r,,(g), T(g-i, g-2) = S,,5kiT,jkliqi,q2), (A.2) 
T(^)(g-i,g-2) = ^5.iei?(-|2)T,,fcKg-i,92), T^''''HquQ2) = -/:^\-qi)et'\-q2)T,jui{qiA2). 
The various contractions of helicity tensors appearing in the main text are 

eS^«)(g-) = vr,,(g-i)6lf )(-!,), = ^.,(9-2)6;f Vis), 

e(^^^^^3)(^- ) = ^^^^^^^^i)^^^ (A.3) 



where the projection operator vTjj is given in (j2.5p . 

We may explicitly evaluate these contractions in terms of the magnitudes qi of the momenta 
and the helicities Sj by introducing a basis for the helicity tensors. To do so, we first observe 
that the momenta qi lie in a single plane due to momentum conservation. Taking this plane 
to be the (x, z) plane, we may then write 

qi = qi{sm9i, 0, cos6'j) (A.4) 

where the magnitudes qi > 0, and without loss of generality we may choose 9i = 0, < 62 < 
and TT < 6*3 < 27r so that 

n (gj -qI- Q2) ■ n A (gj - qj - ql) . n A 
cos 62 = — , sm 62 = TT^^ , cos 6^3 = — ^ — , sm 6^3 — 



2gig2 ' 2^152' 2gig3 ' 251^3' 

(A.5) 
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with A as given in (|4.2p . The required helicity tensors then follow by rotation in the {x,z) 
plane: 



1 

7^ 



cos 



-iSi sm Ui 
sin^ 6j 



iSi cos t/j 
isi COS 6i —1 
Y — sin 6*4 cos 0j — isjsin^j 

The contractions of helicity tensors used in this paper are then 

„r+i._. A2 



sin 9i cos Oi \ 



(A.6) 



16^/2c 



e(++-)(g-.) 



4\/262 



23 



A2 



123 



16^/2( 



-fe - 012)' 



123 



^(++)(- ) ^ «123(«23 -^-l) 



8&i3 



(+-) 



in 



(«13 - g2)^(Ql2 - 93) 

8bh 



^(++),_^ 0123(023-91) rr,-2- /- -N \21 

= T^^2 [2^10123(023 - gi) - A J, 



16c 



123 



(oi3 - 92)(oi2 - 93) r^-2/- - w- -n , 721 

^(^i) = [2gi(ai3 - g2)(ai2 - 93) + A J 



16cf23 



(A.7) 



Let us now discuss the spinor helicity formalism introduced in [7] . For any three-dimensional 
vector q we can consider a four-dimensional vector q'^ = (q, q) which satisfies q^q^ = 0. There- 
fore, q can be represented by spinors A^ as q^ = c^^A^A", where we use the same conventions 
as in [7], namely 



1 
-1 



q'^q2^. = -2(12)(l2), 



(12) — AiaAf — eafeA^Ag, (12) — AidA2 — e^j^A^Aa. 



(Ai 



Here, a'^"^ is a vector of Pauli matrices. The spinors Aia and A2a corresponding to the two 
momenta qi and 52 are denoted by |1) and |2), respectively. Spinor indices are raised and 
lowered by means of e"^^ and its inverse. 

To compare our results with those of [7j, we need an explicit expression of the inner 
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products (12), etc., in terms of momenta. A possible solution for a spinoio is 



where the second expression makes use of the spherical coordinates 

q = g(sin6'cos0, sin^sin^, cos^). (A. 10) 

Choosing momenta qi and q2 as in ()A.4p . and making use of ()A.5p . we find 

(12) = -V^sin =-\yJ{qi+q2Y-ql (A.ll) 

In general the sign depends on the orientation of (^1,^2)- Since we choose ^ < 62 < tt, the 
orientation is assumed to be positive. Note in particular that (21) = —(12). Similarly, we find 



1 \ 77 ,— , 1 



(12) = -^ql - {q, - q2Y , (12) = --^{q,+ 92)2 - eg = (12) , (A.12) 
in agreement with (B.6) of [7J. Combining these results we find 

[(12)(23)(31)]2 = ^a?23 = -^cl23e(+++)(g-.), 

[(12)(23)(31)]2 = ^(ai2-g-3)' = -^430^++"^fe), (A.13) 

which we have made use of in the main text. 

In four dimensions, the complexified symmetry group is locally isomorphic to SL{2, C) x 
5L(2,C), in which case dotted and undotted indices transform independently. In our case, 
however, the symmetry group is that of rotations of three-dimensional space, which corre- 
sponds to SL{2,C) ^ SL{2,C) X 5L(2,C) embedded diagonally. An additional invariant 
tensor therefore exists, which we may choose to be e"°: this means that we are now allowed 
to contract dotted with undotted indices. In particular, 

(AA) = WaaX"" = -X^iX")* = -X\~Xy = -q, (A.14) 

motivating the following definition for complex conjugates 

Aa = (A*^)* = -e,aA^ A"^ = {XaT = -e'^'^A,. (A.15) 

The helicity tensors used in [7j may now be defined as 

{s)abab _ Js)aaJs)bb /' A 1 



^We corrected signs in (B.l) in v2 of [7] so that g'' = cr^^A^A" 
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where 



/.(+)ad 



A" A" 



A-)aa 



Contracting with Pauh matrices, we then find 
/ \ 



COS u COS (p — I sm ( 
cos 9 sin (j) -\- i cos < 
— sin 6* 



A'^A'^ 



(A.17) 



'MP ~ ^ 







COS 6 cos + i sm ( 
cos 9 sin (/> — i cos ( 
— sm9 



(A.18) 



Since the time components vanish, we may regard the C^^jT as three-dimensional vectors, and 
if ^/i^M^f = in four dimensions, then clearly ^i^^p = in three dimensions as well. 

Let us now compare the vectors Cmj^ with those implicit in our own convention ()A.6P for 
the helicity tensors. In our case, we started with 

= (0,l,zs,0) for = (gi,0,0,gi), (A.19) 

and then obtained all other ^^^'^^{q) by rotation in the {xz) plane. In this way, we find that 



?MP V 



(A.20) 



Our normalisation of helicity tensors is then such that 



= -^^MPiii-^- 



(A.21) 
(A.22) 



thus we find 

e^'Va) = 

^/2 ^^^^ '''' '' 

Finally, we chose to define T('\q) = \e^'^ (-^Ti.j{q) (so that Tij{q) = T^'\q)e\f {q)), which 
leads to 



1 



2V2 



(A.23) 



where [7| defines instead T^p{q) = e^p {q)Tij (q) . 



B Evaluation of integrals 

To evaluate the holographic formulae (|3.4p - (j3.6p we must compute specific helicity projections 
of the stress tensor 3-point function. One option, discussed in Section [B. II below, is to project 
into a helicity basis at the very outset of the calculation, leaving only relatively straightforward 
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scalar integrals to evaluate. A second option, discussed in Sections IB. 21 and IB.31 is to directly 
evaluate the tensor integrals for the full 3-point function, projecting into a helicity basis only 
as the final step. While the evaluation of tensor integrals is more demanding, the Ward 
identities permit a useful consistency check of the results. The required evaluation of tensor 
integrals may be accomplished using either a method due to Davydychev |32t [33], or else via 
a Feynman parametrisation approach, as discussed in Section IB. 21 and IB. 31 respectively. In 
practice, we computed integrals using all three methods and cross-checked the results for each 
method against those of the others. 

B.l Helicity projection to scalar integrals 

To illustrate the steps involved, let us consider the following integral derived from the result 
(j4.7p for minimal scalars 

^ J q^iq + qiriq-qsr 

Making use of the explicit basis ()A.6P , we find 

V2elf\-q3)qiqj = ql cos^ Oz+fz G^i-q^-lq^qz sin 6*3 cos e3+2is3qyqz sin 6'3-2is3g^gj, cos 6*3. 

(B.2) 

Since the external vectors iji all lie in the {x,z) plane and thus have no y-component, the 
imaginary part of the integral (]B.ip is odd under qy — )• —ijy and therefore vanishes. To deal 
with the remainder, it is then convenient to replace q^ = (f' — q^ — q^ and to substitute for q^ 
and qx according to 

q- qi _ 1 _ _ cot ^3 _ _ 
9z = ^— , qx= _ . . g • 93 = — q-qi- (B.3) 

qi qs sm 63 qi 

Here, trigonometric expressions involving 63 are equivalent to specific combinations of external 
momenta according to ()A.5p . Finally, using the standard replacements 2q ■ qi = (g + gi)^ — 
q^ — qf, etc., the integral (jB.ip may be reduced to a sum of elementary 2-point integrals and 
a single 3-point integral, 

^"^"^k^iq + qiyiq-QsV^ml^i ^^'^^ 
(Note that this latter integral reduces to a standard 2-point integral upon substituting q' = 

q/f and q[ = qi/qf). 

The evaluation of all remaining helicity-projected 3-point integrals proceeds in a similar 
fashion, the only complexity arising from the need to keep track of moderately large expres- 
sions. 
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B.2 Tensor integrals via Davydychev recursion 

An elegant general method for evaluating tensor Feynman integrals corresponding to arbi- 
trary 1-loop A^-point diagrams was proposed by Davydychev in |32^ 133]. Here, we review its 
application to the tensor integrals appearing in our calculations of the stress tensor 3-point 
function. 

Our goal will therefore be to evaluate massless 1-loop 3-point integrals of the general form 

where, for reasons that will be apparent shortly, we have kept the spacetime dimension n, as 
well as the powers (where i = 1 ... 3) appearing in the denominator, arbitrary. We will 
temporarily denote spacetime indices with Greek letters to avoid confusion with the index i. 
Note also that our choice of Euclidean signature will result in a few minor changes with respect 
to the corresponding formulae reported in jSJUHsjil. The symmetric form of the momenta in 
the denominator is convenient; to recover the form of the momenta in the denominator used 
in the main text one simply shifts q^.^ {q — P3)fi (see Fig. H}. 

In [32], a general formula was derived allowing the tensor integral (jB.SP to be expressed as 
a sum of symmetric tensors constructed from the spacetime metric and the external momenta, 
multiplied by coefficients given in terms of scalar integrals of the form 

J{n-{vi\)= I dJ'q- -T5— ^ , -75—. (B.6) 

Explicitly, in Euclidean signature, this formula reads 



^ — ' Z vr 



X {[g]^[pir[P2rHP3r}^,...^.MJin + 2(M - A); {z., + K,}), (B.7) 

where (z/)^ = r{v -\- K)/r{v) is the Pochhammer symbol and [^2]''^ [Psl'^^l/ti-.-MA/ 

denotes the symmetric tensor constructed out of A copies of the metric tensor and Kj copies of 
each momenta pi. (Thus, for example, {fi'Pi}Mi/'2M3 = 9^ll^l2Pl^l^+ 9^ll^l^Pl^l2+ 9^l2^l^Pl^ll^ where 
for present purposes the metric tensor g^y = 5^^ ) In the formula (IB.Tp . the sum runs over all 
possible non- negative values of A and Kj, such that the total rank 2A + equals M. Note 

in particular that the values of n and Vi appearing in the scalar coefficient integrals differ from 
those appearing in the original tensor integral ()B.5p . 



^^Alternatively, one could use the Lorentzian formulae quoted in [321133) and continue to Euclidean signature 
after completing all computations. 
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Figure 4: Labelling of momenta 



Equipped with the general formula ()B.7p . we may therefore reduce tensor integrals of the 
form (IB.SP to scalar integrals of the form (IB.6j) . The evaluation scheme is then completed 
by a set of recursion relations enabling the scalar integrals ()B.6P to be reduced to elementary 
integrals. In [33j, it was shown that 



J(n; {i^i,U2, J^3 + 1}) = ^ ^2-2 k(2z^i + z^2 + 1^3 - ri)qj 

+ {2u2 + J^i + 1^3 - n)ql - {2u3 + ui + 1^2 - n)ql)j{n; {1^1,1^2, J^s}) 

+ i'2qlJ{n\{h'i -1,1^2 + 1,1^3}) +viqlJ{n;{vi + l,V2- l,z^3}) 

+ v^qlJin; {ui - 1,U2, 1^3 + 1}) - J^iqi-Ji^'^ {^^i + 1^ ^^2,2^3 - 1}) 

+ i/3g|j(n;{z^i,i/2 - l,z^3 + 1}) - z^2^i'/(n;{i^i,i^2 + l,z^3 - 1}) , (B.8) 

with similar formulae for J(n; z^2i 1^3}) and J(n; {1^1, i^2+li 1^3}) following by permutation 

of indices. If we regard the indices (1^1, i^2) ^^3) as coordinates on an integer lattice, these 
recursion relations allow us to construct three integrals in the plane Yli i^i = c + 1 in terms 
of six contiguous integrals in the plane Yli — ^- Now, in general, we are interested in the 
region Vi > 0. Any integrals on the boundary of this region may be evaluated trivially: if 
more than one of the vanish the integral is zero in dimensional regularisation, and if only 
one of the Vi vanishes, the integral reduces to the standard 2-point integral 

T{iyi)T{u2)T{n - vi - V2) 

The first non-trivial integral for which all the ui > is therefore J(n; {1, 1, 1}), which sits in 
the plane cr = 3. From this integral, plus the appropriate 'boundary' integrals, we may then 
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use the recursion relations (jB.SP to construct all the non-trivial integrals in the plane a = A, 
namely J(n; {1,1, 2}) and its permutations. Through repeated appHcation of the recursion 
relations, we may proceed to evaluate any integral with positive integer in terms of the 
initial integral J{n; {1, 1, 1}) and boundary integrals of the form (jB.9p . 

Examining the form of ()B.7p . we see that to evaluate a tensor integral of rank M in 
three dimensions, we need to evaluate the corresponding scalar integrals J{n; {z/j}), and hence 
initial integrals J(n; {1, 1, 1}), in all odd dimensions 3 < n < 2> + 2M. As noted in the previous 
subsection, the initial integral in three dimensions may be evaluated by inverting all momenta 
yielding 

J(3;{1,1,1}) = ^^. (B.IO) 

919293 

The higher odd-dimensional initial integrals may then be obtained using the additional recur- 
sion relation 



J(n + 2; {1,1,1}) 



27r 



W2+ir(n/2) 



nvr 



COSeC — {Hn-2H2 - 2Hn) 



{n-2){2Hi- Hl)l r(n-l) V2 

+ (9iM3)V(n; {1,1,1})], (B.ll) 

where 

Hn^q^ + q^+q^, (B.12) 

which may be derived from (jB.Sp and the Euclidean analogue of equation (6) in |32| . 

Armed with the above analysis, the evaluation of tensor integrals of the form (IB.Sh is now 
straightforward. For the computation of stress tensor 3-point functions in the main text, we 
need to evaluate the six tensor integrals with ui = 1^2 = 1^3 = ^, n = 3 and ranks M = 1 . . . 6. 
For the lower ranks the calculation may easily be executed by hand, yielding for example 

/ [d9]^27^ — I - ^2 = o- ^- — (^29iM - 9192/.), 

J q [q - qi) [q + q2r 8ai23Ci23 
q/iqu _ 1 



/ 



{ai23Ci23^fj.u + 92(013 + 2g2)9lAt92i. 



9^(9-91)^(9 + 92)^ 16af23Cl23 

+ 9i(«23 + 2qi)q2fj,q2u - &1291m92i/) , (B.13) 

but for the higher ranks it is convenient to automate the process. In the highest rank case 
M = 6, we see from (IB.7P we need to evaluate scalar integrals in planes up to a = 9 (thus 
requiring up to six iterations of the recursion relation (jB.SP ). for odd spacetime dimensions 
up to n = 15. Having explicitly computed all tensor integrals up to rank six, any 3-point 
function {{Tij[qi)Tki{q2)Tmn{q3))) may now be directly evaluated. After checking against the 
Ward identities for consistency, the result may then be projected into the helicity basis. 
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B.3 Tensor integrals via Feynman parametrisation 

In this third method, the correlation functions were again calculated directly in the tensor 
representation ((Tj^jjTjjjj^^sjs)) ^-^d then projected into the helicity basis or traced. All 3- 
point functions we consider may be expressed as a sum of the integrals of the form 

where ti^j^i^j^i^j^ is a tensor build up with q,qi,q2 and a metric 5. In order to calculate this 
integral, Feynman parameters xi,X2,X3, such that xi + X2 + = 1, may be introduced. This 
leads to the substitution q = 1 + X2qi — Xiq2 and the integral takes the form 

[<iq]-^iimmiMi_^ = 2 [ dx [ [dr]^iip22i^, (b.is) 



qHq-qiyiq + q2r J[o,i]^ J ' ' (z^ + a) 



where 



dX = dxidx2dx3(5(j;i + X2 + X3 — 1), 
A = qlx2{l - X2) + q2Xi{l - xi) + 2{qi ■ q2)xiX2 

= qlx2X3 + q2XiX3 + qlxiX2, (B.16) 

and the integration is over the cube {xi,X2,xs) G [0, 1]'^. Finally, we decompose the integral 
(jB.lSp into a linear combination of integrals of the form 

2 1 dXP{x,,X2,x,) I [d/ f"|^-^-^-i^- , (B.17) 

where 2:2, 2:3) is some polynomial in Feynman parameters. The integral over momenta 

may be evaluated by means of the formula 

O f ■ ■ ■ ^irnhm r(3/2 — Ul) S i-^^j-^ \ m-Zl2 /T3 1 g\ 

V [^'1 (P + A)3 = ' (^-'^^ 

where Si^j^,„i^j^ is a completely symmetric tensor constructed from metric tensors with all 
coefficients equal to one. Due to the / i— >• — / symmetry the integrals with an odd number of 
momenta / vanish. 

The remaining task is to evaluate the integrals 

y dXP(xi,X2,X3)A'"-3/2 (B.19) 

over the Feynman parameters. For d = 3, the r.h.s. of (jB.lSP is a well-defined expression for 
any integer m, and ()B.19P exists for any polynomial P and any non-negative m. It turns out 
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that in order to find all the integrals we need of the form ()B.19p , it is enough to evaluate only 
one integral. This integral, coming from six Vs in the numerator of (jB.lTp . is 



/ dXA^/^ = [30^23 - 9ai23&123 + 3a?23^123 + 3a?23Cl23 + 3ai23fel23Cl23 + 2c? 

J D4U a 123 



"l23j 

'123 

(B.20) 



The remaining integrals we need may now be evaluated by the following tricks: 

• Differentiating an integral with respect to ^3 introduces Feynman parameters X1X2, e.g. 



dXxiX2Ai/2 = ^ . A / dXA3/2. (B_21) 



I ^"'^ 3^3 dqs 

Notice that this operation decreases the power of A by 1. 

Integrals such as J dXxf A^/^ cannot be obtained by the above method. In this case, 
we may use the following formulae 

k{ql - ql) I dXxr^A'^-l 

= {qf - <if)B{n + k + l,k + l)- kql J dXx^ (X3 - X2)A''-\ 

= qfB{n + k + l,m + k + l)-kql j dXx>^(x3 - X2)A^-^ 

-mj dXx^x™"^A^ for m > 0, (B.22) 



with numbers k, m and n such that these expressions exist, and where B is Euler's beta 
function. For example, taking k = 3/2, m = 0, n = 1 we find 

93-92 



j dXx?AV2 = _^iL^ j dXxiX2AV2 _ j dXxiX3AV2 



_^ vr g| + 9293 + 93 



128 q2 + 53 

(B.23) 



Integrals with odd numbers of Feynman parameters may be obtained from the integrals 
with even numbers of Feynman parameters by utilising the fact that xi + X2 + X3 = 1. 
For example, 

j dXxi^^''^ = j dXxl^^/^ + j dXxiX2Ai/2^ j dXxix^A^''^, (B.24) 
where the integrals on the r.h.s. may be found in previous points. 
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Iterating the trick described above we may find integrals with different powers of A, e.g. 

I dXA^/^ = Z2 [«123 - ai23&123 - C123] , 



/ 



dXA-V2 = J^, f dXA-3/2 = (B.25) 

ai23 J C123 



This method ahows the exact tensor representation of any 3-point function we consider in 
this paper to be calculated. However, since we are interested in the helicity representation, 
it can be significantly simplified. Any rank-six tensor iiijii2j2j3i3 built out of the metric and 
two independent momenta may be represented as a sum of 499 simple tensors. In general, 
the independent momenta may be different for different tensor indices: we choose to use 
momenta qi and q2 for ii and ji indices, q2 and ^3 for ^2 and j2 and 53,^1 for 13, js- Due to 
various symmetries and Ward identities on t, however, the number of independent tensors is 
usually much smaller. If we consider a 3-point function ((Tajfe^Tajfea^asfea)) projected onto the 
transverse-traceless part, we find only five independent coefficients, i.e., 

^iljiaibi^i2j2a2b2^i3j3ii3b3 ((^£1161-^(1262-^0363)) (B.26) 
= Il^,j,a,b^Ili2j2a2b2^i3j3a3b3 [Ai{qi,q2,q3)S^''"6''^''6''''' 
+ A2{qi,q2,q3)S''''"6'^''-^if,'if2' + (1 ^ 3) + (2 O 3) 
+ A3{qi,q2,q3)5''"''5'''''(t^'q'l' + (1 O 3) + (2 O 3) 
+ A,{q,,q2,q3)S^^'''4'Q3fl'^l + (1 O 3) + (2 O 3) 
+ A,{q,,q2, g-3)9T9t ^T^t 

The coefficients Aj may be easily expressed in terms of the coefficients of various tensors 
appearing in {{Ti^j^Ti^j^Ti^j^)) . Specifically, we see that 

^1(^1,92,93) = 8- coefficient of (5^i^2 5^2i3 5i3ii^ 
^2(91,92,93) = 8- coeflicient of <5*i^2 5^2i3^i3^ii^ 
^3(91,92,93) = 2- coefficient of <5^i ^2 5^ii2g»3^J3^ 

^4(91,92,93) = 4- coefficient of <5^i ^2 g^igfgi'9f, 
^5(91,92,93) = coefficient of 92' 9? 9^9? 9f- 

In other words, it is enough to calculate only five scalar integrals in order to evaluate the five 
independent coefficients Aj. 

Finally, to obtain the result in the helicity basis we may contract ()B.26P with helicity 
tensors. Using the identities (jA.ip and ()A.3P one finds five contractions of helicity tensors 



44 



corresponding to the independent transverse-traceless tensors. These results were checked 
by a simple computer algebra program which carried out a brute force calculation of all 
499 coefficients in ((Ti^jjTjjjjTjgjg)) before projecting the result into transverse-traceless and 
helicity bases. (This procedure also enables checking against the Ward identities.) 

Note this method also works if some indices are traced. In this case, the situation is 
analogous to that for tensors of ranks two and four. Transverse-traceless tensors of rank four 
have three independent coefficients, while those of rank two have only one. The coefficient 
can be evaluated by the same method as described above. 



C Contribution of ghosts and gauge-fixing terms 

To evaluate the gauge field contribution to 3-point functions we must gauge-fix and introduce 
ghost fields. This procedure generates a new contribution to the stress tensor that depends 
on the gauge-fixing part of the Lagrangian. Here we show that this part does not contribute 
to the 3-point functions. The general argument is based on the fact that the full Lagrangian 
for the gauge field is 

Sym = -^ I d^x tr \-MfI^ + 5bO\ , (C.l) 

where O is a gauge-fixing part containing ghosts and 5b is an infinitesimal BRST transforma- 
tion. The full stress tensor is therefore 

T^^ = T,^ + Tfl (C.2) 

where is a BRST-exact operator. Since physical states correspond to the cohomology of 
the BRST transformation, Tf^ vanishes when acting on such states. Therefore, inside any 
vacuum correlation function, can be replaced by T^. 

As this is a formal argument, we will also present now an explicit perturbative proof that 
the gauge- fixing part does not contribute to any correlation functions. We work in the 
gauge and to first order in ^ym- '^^^ ghost part and gauge-fixing part of the action may be 
written as 



9ym 



[ d^x tr l^B^f + A^diB^] , 5gh = [ d^x tr [diC^diC^] . (C.3) 

where and are the antighost and the ghost fields, and is the BRST auxiliary field. 
All fields are in the adjoint representation and are regarded as traceless hermitian matrices. 
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The full Yang- Mills theory is given by the action 



Sym 



1 



9ym 



d^x tr 



1 



?I J7l 



F 

4 u «J 



+ S"^ + S'gh- 



This leads to the following propagators 

{{B\q)B\-q))) = 0, 

and 



iiB\q)F^'A-m = 0, 



- {{Al''{q)B^\-q))) = {{{diC^''){q) c-^\-q))) = J'^*^^-^ '^YM^i 



r 



(C.4) 

(C.5) 
(C.6) 



Here, by {diC^"'){q), we denote the Fourier transform of diC^°'{x). 

The stress tensor and the T tensor defined in (j2.2p corresponding to each component of 
the action is given by 

1 .„r_r . 1 



^ ijkl 



2 ^^[Ki^jk - ^ijjHlHl]^ 

9ym ^ 

tr[-P,jkiAidiB' + (^i4(^')']' 
9ym ^ 



1 



5ym 



^ tr[-6,^kSi)jAl^d^B' + 6ijAl^di)B' - 6,^kSi)A{B'f]d{x - y), 
9ym ^ 



^ ijkl 



9ym 

where Qijkimn is defined in (j4.36p . The full stress tensor and T tensor is a sum 



rpYM _ rpA I rp^ , rpgh ^YM _ , 'yS. , 'V^gh 

^ij ~ -'-ij -^jj ~^ ^ij ' ^ijkl — ^ijkl ^ijkl "i" ^ ijkt 



(C.7) 



(C. 



The mechanism for cancellation of ghost and gauge-fixing terms is very general. Let us 
consider a set of general gauge-invariant operators of arbitrary tensor structure, indexed 
by a, quadratic in field strengths F^ . Consider moreover gauge dependent terms B^"'^ and 
ghost terms C^"^ of the schematic form 

1 



9ym 



9ym 



tr 



tr 



(C.9) 
(C.IO) 



where of'^"^ is linear in B^ , Q^^i'^) ig quadratic in B^ and O^'^"^ is linear in , but are 
otherwise operators of arbitrary tensor structure which may contain derivatives but no other 
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fields. We consider operators O^"^ = + B^"'^ + C^"-* and their n-point function in the 
Yang- Mills theory with the action ()C.4p . The stress tensor and the T tensor are of this form. 
We find 

+ (0(1)^(2) ^ ^ ^ ^(n)^ ^ (^(1)0(2) . . . ^W) + . . . + (^(1)^(2) . . . gW) 
+ . . . 

+ ^ ^ ^ ^(n)^ ^ (C(1)C(2) . . . C(")) (C.ll) 

since there is no interaction between ghosts and any other fields at leading order in 5ym- 
will now show that all terms but the first one cancel. 

To begin, we observe that all terms containing at least one and at least one B vanish. 
Indeed, when Wick's theorem is applied, there must be at least one contraction between F 
and B fields, or between B and another B field, which gives zero by (jC.Sp . 

Now consider the term with B operators only. When expanded, it has 2" terms, but every 
term containing {B^)^ must evaluate to zero as there must be at least one B-B contraction. 
Therefore, only one term survives, namely 

(^(1) . . .^H) = ^(tr (40//^)[i?^^]) • . . . • tr (4::6/j(")[i?^"])). (C.12) 
9ym ^ ^ ^ ' 

The only non-vanishing way of contracting fields is to contract auxiliary fields with gauge 

fields. This gives precisely the same possible contractions as in the ghost part, which is 

. . .CW) = 4r(t^ feiC-^^Oj(^)[c^^]) • . . . • tr fe„c-^"0£(")[c^"])). (C.13) 

It follows that if -Of '^"^ = Of'^"^ for ah q, then ([UT2]) and ([UT3l) cancel each other out, 
due to (IC.6P and the anti-commuting nature of ghost fields. 

In our case, and to this order in (7yM' there are no gauge boson interactions and so the 
gauge group G is effectively ?7(1)^™^. For (i;^^ 7;YM j^yM) therefore find 

= T^. , B'' = T^ . , = , (C.14) 

and 
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for a = 1, 2, 3. For (T^.^ T^^), we have 



ijkl 



h(2) _ 

— -'■mi 



and 



AA,(1) 
ijkl,m ijkl,m 



-o. 



A,(2) ^ AC,(2) ^ 
van.i mn.i 



Pmnik^k- 



(C.16) 



(C.17) 



It follows that the contribution due to the gauge-fixing part of the action indeed cancels out. 



D Further Ward identities 

D.l Conformal Ward identities 

The conformal Ward identities are given by 







{i ^ j) 



((Oi(gi)02fe)03fe))), (^,i=l,2,3) (D.l) 



2(A2 - d)d2^ + {-npA + 



((r,,,,(g-i)02fe)03fe))) 



+ 2 [iSi,i,d^' - + (ai o h,ii o ii)] ((r„,fe,fe)r,,,, (92)0(^3))) 

+ 2 [(J,,^^^"^ - <J^^52,,)55^^ + (a2 O 62,^2 ^ ^2)] m^,M)Ta,bM)0{q3))) , 
2 



{{Thji iQi)Ti2j2 {Q2)Ti^n (^a))) 



+ 2 [(<5ii^a?i - + (ai o 61, zi o ji)] {{TaMQi)TnniQ2)Ti,jM))) 



+ 2 



- '^^^52^J5;^^ + (a2 ^ 62,Z2 ^ ^2) ((T,,,, (ft )r„,b, (g-2)T,3,3 (g-3))) , 



(D.2) 



(D.3) 



(D.4) 



where Oj is taken to have dimension Aj and = d/dqf. 
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D.2 DifFeomorphism Ward identity 

The diffeomorphism Ward identity for 3-point functions may be evaluated by functionally 
differentiating \7^{Tij{x))s = twice with respect to the metric pi]) yielding 

= qii{{Tij{qi)Tkl{q2)Tmniq3))) - '^qii{(^ijmn[qiA'i)Tkl{q2))) - 2gii((Tjjfc;(gi 

- '2qi(k{{Tl)j{q3)T.^n{-q3))) - '2qi(m{{Tn)j{q2)Tkl{-q2))) - Skiq2p{{Tpj{q3)Tmn{-q3))) 

- Smnq3p{{Tpj{q2)Tki{-q2))) + q2j{{Tki{q3)Tmn{-q3))) + q3j{{Tmn{q2)Tki{-q2)))- (D-5) 

We explicit checked that all our 3-point functions satisfy this identity. Note that our result 
differs from that quoted in |14] due to a difference in the definition of the 3-point function: 
here, we define the 3-point function by the insertion of three copies of the operator T,j, whereas 
in [H], the 3-point function is defined via functionally differentiating the generating functional 
three times. These two definitions differ from each other by semi-local terms (see the discussion 
around (f3TT]) ) . 
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